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ABSTRACT 

O ' We use the catalogue of galaxy groups constructed from the Sloan Digital Sky Survey 

(SDSS DR4) by Yang et al. (2007) to study the dependence of galaxy activity on 
^ ' stellar mass, M*, halo mass, Mh and group hierarchy (central vs satellite galaxies). 

The wealth of data provided by the SDSS allows us to split the sample on the basis of 
galaxy activity in star- forming galaxies, galaxies with optical AGN activity, composite 
galaxies (both star formation and optical AGN activity) and radio sources. We find a 
^ \ smooth transition in halo mass as the activity of central galaxies changes from star for- 

^P . mation to optical AGN activity to radio emission. Star-forming centrals preferentially 

[^ I ' reside in haloes with Mh < 10^^ /i~^M0, central galaxies with optical- AGN activity 

typically inhabit haloes with Mh ^ 10^'^ /i~^Mq, and centrals emitting in the radio 
mainly reside in haloes more massive than lO^'*ft.~^M0. Although this seems to suggest 
r~^ ' that the environment (halo mass) determines the type of activity of its central galaxy, 

^^ . we find a similar trend with stellar mass: central star formers typically have stellar 

masses less than \Q^^h~'^ M©, while optical- AGN hosts and central radio sources have 
characteristic stellar masses of ^ 10^"'^/i~^ M© and ^ lO^^'^ft."^ Mq, respectively. 
Since more massive haloes typically host more massive centrals, it is unclear whether 
the activity of a central galaxy is causally connected to its stellar mass or to its halo 
^ , mass. In general, satellite galaxies have their activity suppressed with respect to cen- 

jrt ■ tral galaxies of the same stellar mass. This holds not only for star formation activity, 

but also for AGN activity in the optical and the radio. At fixed stellar mass, we find 
that the activity of satellite galaxies depends only weakly on halo mass. In fact, using 
a set of reduced conditional probability functions, we find that for satellite galaxies 
the dependence of galaxy activity on halo mass is more than four times weaker than 
the dependence on stellar mass. We also investigate whether the strength of galaxy 
activity is infiuenced by the environment, using the specific luminosities in Ha, [OIII] 
and at 1.4 GHz. We find that the specific luminosity of star-forming and optical- AGN 
centrals systematically decreases with increasing halo mass, while the specific lumi- 
nosity of centrals with high radio activity increases with halo mass. Independently 
of the activity class, the specific luminosity of satellite galaxies does not significantly 
vary with halo mass. All these results are consistent with a picture in which low mass 
haloes accrete cold gas, while massive haloes have coronae of hot gas that promote 
radio activity of their central galaxies. 

Key words: galaxies: active - galaxies: clusters: general - cosmology: dark matter - 
radio fines: galaxies. 
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1 INTRODUCTION 

Ever since the seminal work of Hubble (1926), it is known 
that the local galaxy population consists of two main 
branches: the red sequence populated by bulge-dominated 
galaxies with old stellar populations (early- type galaxies), 
and the blue cloud comprising disk galaxies with some level 
of star- formation activity (late-type galaxies). With the help 
of large galaxy redshift surveys, such as the Sloan Digital 
Sky Survey (SDSS; York et al. 2000; Stoughton et al. 2002), 
this bimodality has been dramatically confirmed (e.g., Blan- 
ton et al. 2003; Strateva et al. 2001). This bimodality is 
tightly linked to environment, with early-types preferen- 
tially residing in dense environments such as galaxy clusters, 
while late- types mainly populate low-density structures such 
as galaxy groups and the more generic field (e.g.. Dressier 
1980; Postman & Geller 1984; Lewis et al. 2002; Goto et 
al. 2003; Kauffmann et al. 2004; Hogg et al. 2004; Wein- 
mann et al. 2006). 

In the current paradigm of galaxy formation, it is be- 
lieved that virtually all galaxies initially form as late-type 
disk galaxies due to the cooling of gas with non-zero angular 
momentum in virialized dark matter haloes (e.g.. White & 
Frenk 1991). During their subsequent evolution, disk galax- 
ies may be transformed into early-types via a variety of 
transformation mechanisms. Probably the most important 
of these is major mergers, in which two disk galaxies of 
roughly equal mass merge to produce a spheroidal galaxy 
(Toomre & Toomre 1972; Negroponte & White 1983). Dur- 
ing such a merger the gas looses its angular momentum and 
accumulates at the center where it can trigger an intense 
starburst and fuel the central black hole (e.g., Hernquist 
1989; Mihos & Hernquist 1996). This not only couples the 
growth of the central black hole to that of its host galaxy, 
which may explain the observed Mbh - o" relation (Geb- 
hardt et al. 2000; Ferrarese & Merritt 2000), but it may 
also quench subsequent star formation by exhausting the 
gas supply and by expelling or heating the cold gas via 
AGN activity (e.g., Menci et al. 2005; Springel, Di Matteo 
& Hernquist 2005; Hopkins et al. 2006). In addition to this 
so-called quasar mode of AGN feedback, one also considers 
the so-called radio mode, in which the radio activity associ- 
ated with a low accretion efficiency of the central black hole 
provides a heating term that can offset the cooling, there- 
with preventing any new gas from cooling and forming stars. 
It has been suggested that this radio-mode feedback plays 
a crucial role in explaining why the most massive galaxies 
are red and dead (e.g.. Bower et al. 2006; Croton et al. 2006; 
Cattaneo et al. 2006; Kang, Jing & Silk 2006; Nusser, Silk & 
Babul 2006; Sijacki & Springel 2006). The inclusion of AGN 
feedback (radio mode and/or quasar mode) allows models 
of galaxy formation to reproduce the bright end of the ob- 
served luminosity function. In order to reproduce the faint 
end of the luminosity function, the models typically invoke 
supernova feedback to expel gas from low mass haloes (e.g., 
Larson 1974; White & Rees 1978; Dekel & Silk 1986). 

Clearly, star formation and AGN activity are thought to 
play a pivotal role in the formation and evolution of galaxies. 
An important, largely outstanding question, is how this ac- 
tivity is influenced by the galaxy's environment. Such a link 
will put strong constraints on the mechanisms responsible 
for triggering and quenching such activity. For example, if 



AGN were to be fueled by the same cold gas that feeds star 
formation in disk galaxies, their environmental dependence 
would most probably follow the one of star-forming galaxies, 
i.e. the fraction of AGN would decrease in denser environ- 
ments. Alternatively, if the AGN activity were maintained 
by galaxy-galaxy interactions, it would be preferentially as- 
sociated with intermediate mass groups with low velocity 
dispersions or with young clusters that are not yet virialized 
(e.g., Gunn 1979; Shlosman et al. 1990). 

The existence, in the local Universe, of a star-formation 
rate - density relation is well established by now. Hashimoto 
et al. (1998) measured the environment in terms of local den- 
sity, which is based on the distance of a galaxy to its third 
nearest neighbor. They found a continuous correlation be- 
tween star-formation rate (SFR) and local density, whereby 
galaxies with a higher SFR reside preferentially in lower den- 
sity environments. In addition, the authors showed that low 
levels of star formation depend more strongly on the local 
density within a cluster than in the field, while high SFRs 
(which they identified as starbursts) are less sensitive to local 
density in both clusters and the field. This suggest that two 
processes determine the SFR - density relation: gas removal 
via stripping, responsible for the low SFRs of galaxies in 
dense environments, and galaxy-galaxy interactions, causing 
the prevalence of starbursts in intermediate-density environ- 
ments. Similar results were obtained by Carter et al. (2001) 
and Gomez et al. (2003) who found that the fraction of star- 
forming galaxies decreases with increasing galaxy density 
(derived from the distance of a galaxy to its 10th nearest 
neighbor). In particular, Lewis et al. (2002) and Gomez et 
al. found a break in this trend at a critical density of 1 h^^ 
Mpc^, below which the fraction of star- forming galaxies is 
high and constant and above which it drops rapidly. The 
authors also found that the SFR increases with the cluster- 
centric radius, to reach the SFR typically measured in the 
field at radii larger than 3-4 times the cluster virial radius. 
They concluded that this trend can not be explained by the 
morphology - density relation alone, against the intuitive in- 
terpretation of the SFR - density relation as a direct prod- 
uct of the environmental dependence of galaxy morphology. 
Kauffmann et al. (2004) found a decrease in the specific star- 
formation rate (i.e SFR/M.) of a factor of 10 between low- 
and high-density environments. More recently, Weinmann et 
al. (2006), using a catalogue of galaxy groups from the Sloan 
Digital Sky Survey (SDSS DR2), have been able to quan- 
tify the galaxy environment in terms of halo mass, i.e. the 
mass of the dark-matter halo hosting a galaxy group. The 
authors have shown that the fraction of late-type galaxies 
(actively forming stars) smoothly increases with decreasing 
halo mass, down to 10^^ M0. For any given halo mass, this 
same fraction increases with the halo-centric radius. 

As for AGN, a number of studies in the literature have 
found a constant fraction of optical AGN at low redshift 
across a wide range of environments, from the field to galaxy 
clusters (Monaco et al. 1994; Coziol et al. 1998; Shimada et 
al. 2000; Schmitt 2001; Gomez et al. 2003; Miller et al. 2003). 
Kauffmann et al. (2004) obtained rather opposite results, in- 
dicating that the fraction (at fixed stellar mass) of AGN with 
strong [OIII] emission doubles in low-density environments 
with respect to the denser ones. In their study, environment 
is primarily based on the number of neighboring galaxies, 
but a comparison with cosmological simulations allowed the 
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authors to assign to the lowest/highest local densities a char- 
acteristic halo mass of 10^^ - 10^^ Mq and lO" - 10^^ Mq, 
respectively. 

A different approach to studying the environment de- 
pendence of optical AGN activity has been used by Wake 
et al. (2004), Groom et al. (2005) and Li et al. (2006), who 
measured the cross-correlation between optical AGN and 
a reference sample of inactive galaxies. The authors agree 
that optical AGN and the control sample share the same 
clustering amplitude on scales larger than a few Mpc. In the 
range between 100 kpc and 1 Mpc, however, optical AGN 
are clustered more weakly than the control galaxies, while 
for scales smaller than 70 kpc there is a weak indication for 
optical AGN to cluster more strongly than inactive galax- 
ies. Using N-body simulations, Li et al. reveal that the AGN 
anti-bias measured between 100 kpc and 1 Mpc is possibly 
due to AGN residing preferentially at the center of their 
dark- matter haloes. In addition, the fraction of optical AGN 
located in 10^^ - 10^^ Mq haloes is found to be larger than 
that in haloes as massive as lO^'* - 10^^ Mq. 

Based on their [QUI] luminosity, optical AGN are of- 
ten divided in LINERs and Seyferts, where the latter have 
higher [OIII] luminosities and accretion rates (e.g., Kewley 
et al. 2006). Kelm et al. (2004) and Constantin & Vogeley 
(2006) have shown that LINERs are more strongly clustered 
than Seyferts. This trend is consistent with the morphology- 
density relation, given that LINERs are generally hosted by 
more massive galaxies (cf. Kewley et al. 2006). However, 
Constantin et al. (2007) found that LINERs associated with 
moderately massive host galaxies (< 3 x 10^'^ Mq) appear 
to be more common in voids (the most underdense envi- 
ronments in the Universe) than in walls, although their ac- 
cretion rates and old stellar populations are comparable to 
their wall counterparts. Void AGN hosted by more massive 
galaxies seem, instead, to accrete more strongly and to have 
younger stellar populations than those in walls. These trends 
would thus suggest that the AGN activity in voids is tightly 
connected with the availability of fuel, the same fuel that 
supports the star-formation activity of their host galaxies. 

On the basis of their radio activity, AGN are commonly 
distinguished between radio-quiet and radio-loud. These two 
groups are known to exhibit different clustering properties, 
which are also different from those of optical AGN. For 
example, Yee & Green (1984) and EUingson et al. (1991) 
found that radio-loud AGN are located in galaxy clusters 
as rich as Abell class 0/1, possibly reflecting the picture 
where radio-loud AGN reside in elliptical galaxies, while 
radio-quiet AGN are mostly associated with disk galaxies. 
On the contrary, McLure & Dunlop (2001) found that the 
clustering amplitudes of radio-quiet and radio- loud AGN are 
rather similar, and both groups of radio AGN are detected in 
environments as rich as Abell class 0. Later on, Prestage & 
Peacock (1988), Hill & Lilly (1991), Miller et al. (2002) and 
Best (2004) claimed that radio-loud AGN appear to favor 
galaxy groups and weak clusters and to avoid the densest 
environments unless they are hosted by the central cluster 
galaxy. A more in-depth analysis by Best (2004) has shown 
that the fraction of radio-loud AGN with absorption lines 
(in their optical spectra) increases in denser environments 
while the fraction of radio-loud AGN with emission lines 
decreases. Signiflcantly larger samples of radio AGN such 
as those drawn from the SDSS survey have conflrmed the 



early results obtained by Best (2004) and have indicated 
that radio-loud AGN are more likely hosted by the bright- 
est group or cluster galaxies (Best et al. 2005a, 2007; Kauff- 
mann, Heckman & Best 2008). 

In the studies mentioned above, the environment is usu- 
ally parameterized in terms of a projected number density 
of galaxies, above a given magnitude limit. Typically this 
number density, indicated by E„, is measured using the pro- 
jected distance to the nth nearest neighbor, with n typically 
in the range 5-10. However, as discussed in Weinmann et 
al. (2006), the physical interpretation of E„ itself depends 
on environment: in clusters, where the number of galaxies 
is much larger than n, S„ measures a local number density, 
which is a sub- property of the cluster (i.e., E„ is strongly 
correlated with cluster-centric radius). In low-density envi- 
ronments, which are populated by haloes which typically 
contain only one or two galaxies. En, measures a much more 
global density, covering a scale that is much larger than the 
halo in which the galaxy resides. This ambiguous, physi- 
cal meaning of S„ severely complicates a proper interpre- 
tation of the various correlations between environment and 
galaxy properties. Furthermore, this environment indicator 
can not distinguish between central and satellite galaxies. 
This is important, since models of galaxy formation and 
evolution predict that these two types of galaxies interact 
with their host haloes in very different ways: while central 
galaxies can accrete halo gas that cools, and may drive sub- 
sequent star formation and/or AGN activity, satellite galax- 
ies are subjected to various processes such as strangulation, 
ram-pressure stripping and galaxy harassment (Bower 2006; 
Croton et al. 2006; Kang, Jing & Silk 2006; Hopkins et al. 
2006; Larson, Tinsley & Caldwell 1980; Balogh, Navarro & 
Morris 2000; Gunn & Gott 1972; Farouki & Shapiro 1981; 
Moore et al. 1996). 

In this paper, we re-visit the environmental dependence 
of galaxy activity (star formation, optical AGN activity and 
radio emission), using a large galaxy-group catalogue con- 
structed from the SDSS. This catalogue allows us to param- 
eterize the environment in terms of the halo mass in which 
each galaxy resides and to split the galaxy population in 
centrals and satellites, thus allowing for a more straightfor- 
ward comparison with galaxy formation models. In partic- 
ular, the observational results presented here place impor- 
tant constraints on the various mechanisms responsible for 
triggering and quenching different sort of activity in cen- 
tral and satellites galaxies. Throughout this paper we adopt 
a flat ACDM cosmology with Hm = 0.238 and ^a = 0.762 
(Spergel et al. 2007) and we express units that depend on the 
Hubble constant in terms of ft = _ffo/(100 kms~^ Mpc^^). 



2 DATA 

The analysis presented in this paper is based on the SDSS 
DR4 galaxy group catalogue of Yang et al. (2007; here- 
after Y07). This group catalogue is constructed applying 
the halo-based group flnder of Yang et al. (2005a) to the 
New York University Value-Added Galaxy Catalogue (NYU- 
VAGG; see Blanton et al. 2005), which is based on SDSS 
DR4 (Adelman- McCarthy et al. 2006). From this catalogue 
Y07 selected all galaxies in the Main Galaxy Sample with 
an extinction corrected apparent magnitude brighter than 
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r = 18, with redshifts in the range 0.01 ^ z ^ 0.20 and with 
a redshift completeness Cz > 0.7. This sample of galaxies 
is used to construct three group samples: Sample I, which 
only uses the 362356 galaxies with measured redshifts from 
the SDSS, Sample II which also includes 7091 galaxies with 
SDSS photometry but with redshifts taken from alternative 
surveys, and Sample III which includes an additional 38672 
galaxies that lack a redshift due to fiber-collisions, but which 
have been assigned the redshift of their nearest neighbor (cf. 
Zehavi et al. 2002). Since the assignment of optical activ- 
ity to the galaxies (see below) requires SDSS spectra, the 
present analysis is based on Sample I. 

The magnitudes and colours of all galaxies are based 
on the standard SDSS Petrosian technique (Petrosian 1976; 
Strauss et al. 2002) , have been corrected for galactic extinc- 
tion (Schlegel, Finkbeiner & Davis 1998), and have been 
A'-corrected and evolution corrected to z = 0.1, using the 
method described in Blanton et al. (2003). We use the no- 
tation "'^Mr to indicate the resulting absolute magnitudes 
in the r-band. Stellar masses for all galaxies, denoted by 
M«, are computed using the relations between stellar mass- 
to-light ratio and colour of Bell et al. (2003; see Y07 for 
details) . 

As described in Y07, the majority of the groups in our 
catalogue have two estimates of their dark matter halo mass 
Mh: one based on the ranking of its total characteristic 
luminosity, and the other based on the ranking of its to- 
tal characteristic stellar mass. As shown in Y07, both halo 
masses agree very well with each other, with an average scat- 
ter that decreases from ~ 0.1 dex at the low mass end to 
~ 0.05 dex at the massive end. In addition, detailed tests 
with mock galaxy redshift catalogues have demonstrated 
that these group masses are more reliable than those based 
on the velocity dispersion of the group members (Yang et 
al. 2005b; Weinmann et al. 2006; Y07). In this paper we 
adopt the group masses based on the stellar mass rankingj. 
These masses are available for a total of 228062 groups in 
our sample, which host a total of 290729 galaxies that con- 
stitute the sample of galaxies that we use for our analysis 
described below. 

Note that this is not a volume-limited sample. Conse- 
quently, the sample suffers from Malmquist bias, causing an 
artificial increase of the average luminosity (and also stel- 
lar mass) of the galaxies with increasing redshift. To correct 
for this bias, we weight each galaxy by 1/Vmax, where Vmax 
is the comoving volume of the Universe out to a comoving 
distance at which the galaxy would still have made the selec- 
tion criteria of our sample. In what follows all distributions 
are weighted by l/Vmax, unless specifically stated otherwise. 

Using the group catalogue, we split our sample into 
"centrals", which are defined as the most massive group 
members in terms of their stellar mass, and "satellites", 
which are those group members that are not centrals. In 
what follows, we refer to this distinction between centrals 
and satellites as the "group hierarchy". Finally, for each 
group member we determine the projected radius -Rproj from 
the luminosity weighted group center using the angular sep- 
aration and the angular diameter distance at the redshift of 



^ We have verified, though, that none of our results changes sig- 
nificantly if we adopt the luminosity-rank based masses instead. 
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Figure 1. The BPT diagram (Baldwin, Philips & Terlevich 1981) 
for the galaxies in our sample whose [OIII]A5007 and Ha emission 
lines have a S/N ratio ^ 3. The black contours trace 2D bins (of 
0.1 dex xO.l dex) containing an equal number of galaxies (contour 
levels: 3, 50, 500, 1000 and 2000). The dashed line represents the 
demarcation line of pure star formation defined by Kaufl^mann et 
al. (2003b) and the solid line is the extreme starburst demarcation 
line of Kewley et al. (2001). 



the group. We normalize these projected radii by the char- 
acteristic radius of the group, Riso, which is defined as the 
radius inside of which the dark matter halo associated with 
the group has an average overdensity of 180. This radius 
is computed from the mass and redshift of the group using 
equation (5) in Y07. 



2.1 Optical Activity 

In order to assign the galaxies in our sample to one of four 
classes of optical activity, we match the galaxies in our sam- 
ple with the catalogue of emission-line fluxes constructed by 
Kauffmann et al. (2003b). These authors measured various 
nebular fluxes from the SDSS spectra corrected for Galac- 
tic extinction after subtraction of the stellar absorption-line 
spectrum. The latter is fitted with model galaxy spectra 
computed using the stellar population code of Bruzual & 
Chariot (2003), and the best-fitting synthetic spectrum is 
subtracted from the observed spectrum. Emission-line fiuxes 
are then corrected for the intrinsic reddening from the galaxy 
itself using an extinction curve of the form A""'^ (Chariot 
& Fall 2000) and assuming the canonical flux ratio Hq/H/3 
= 2.86 (Case I, Osterbrock 1989). 

We consider a galaxy to be 'optically active' if its 
[OIII]A5007 and Ha emission lines are detected with a 
S/N ratio ^ 3. In order to discriminate between star- 
formation and AGN activity, we use the line-flux ratios 
[OIII]A5007/H/3 and [NII]A6584/Ha to place the galaxies 
in our sample on the BPT diagram (Baldwin, Philips & 
Terlevich 1981). This diagram is traditionally used to sep- 
arate type II AGN from star-forming galaxies and compos- 
ite galaxies. In the standard unifled model (e.g. Antonucci 
1993), AGN are classified as type I if their central black 
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Figure 2. Left-hand panel: the power at 1.4 GHz, Pi. 4, (in W Hz^^) of the radio sources in our sample that reside in groups with an 
assigned halo mass as a function of the absolute magnitude, '^■^Mr — 5Iog(/i), of their host galaxies. The dashed lines indicate the location 
of FRII-type objects. Right-hand panel: Pi, 4 as a function of the galaxy stellar mass, M«. The solid line corresponds to Pi.4/M« = 10^^ 
W Hz~-^/(/i~^Mg ) and separates sources with high radio activity (HRA, P1.4/M, > 10"'^^) from those with low radio activity (LRA, 
Pia/M, < 1012). 



hole and its broad emission-line region are viewed directly. 
When, instead, the central black hole is obscured and only 
its narrow emission-line region is visible, AGN are labeled as 
type II. Since the photometry and spectra of type I AGN are 
dominated by non-thermal emission, they can not be used 
to study the properties of their host galaxies, such as stellar 
mass. Consequently, AGN of type I are excluded from our 
sample. 

Fig. [T] shows the BPT diagram of the galaxies in our 
sample. The dashed line represents the demarcation line of 
pure star formation defined by KaufFmann et al. (2003b, 
their equation 1) and the solid line is the extreme starburst 
demarcation line derived by Kewley et al. (2001, their equa- 
tion 5). We classify as star-forming (SF) all galaxies that 
in the BPT diagram lie below the pure star formation line 
and as type II AGN all galaxies above the extreme-starburst 
line. Those galaxies in between the two lines are labeled as 
composites (CO). The remaining galaxies in sample I that 
either do not satisfy our S/N requirements or do not exhibit 
emission lines are classified as 'not optically active' (NO A). 
The number of star- forming galaxies (SF), optical AGN and 
composite (CO) galaxies is listed in Table 1, where we also 
indicate whether they are centrals or satellites. 

The label AGN comprises both Seyfert II galaxies and 
LINERs (the latter have a lower [OIII]A5007/H/3 at fixed 
[NII]A6584/Ha). Kewley et al. (2006) have demonstrated 
that Seyferts of type II and LINERs form a continuous se- 
quence, where LINERs have a lower accretion rate than 
Seyferts. In ^we will use the specific [OIII] luminosity to 
split our sample in type II Seyferts and LINERS, but for 
the main part of this paper we do not discriminate between 
these two types of AGN. 



2.2 Radio Activity 

In addition to the optical activity, we also consider the ra- 
dio activity of the galaxies in our sample. To that extent 
we match our galaxy sample with the NVSS survey (Na- 
tional Radio Astronomy Observatory Very Large Array Sky 
Survey; Condon et al. 1998) and the FIRST survey (Faint 



Table 1. Overview of Sample 





Centrals 


Satellites 


ALL 


228062 


62667 


SF 


51738 


18017 


CO 


35174 


7358 


AGN 


24273 


4506 


HRA 


2595 


390 


LRA 


2594 


366 



Notes: The numbers of centrals and satellites in our sample and in 
the various classes of optical and radio activity (SF=star forming, 
CO=composite, AGN=optical AGN activity, LRA=low radio ac- 
tivity, HRA=high radio activity). 



Images of the Radio Sky at Twenty centimeters; Becker et 
al. 1995), both conducted at a frequency of 1.4 GHz. We 
follow the identification procedure of Best et al. (2005b), 
and we refer the reader to their paper for details. Contrary 
to Best et al., who only selected galaxies with a flux den- 
sity J5 5 mjy, we select all galaxies detected at 1.4 GHz, 
independent of their flux densities. Using the observed flux 
density and the SDSS redshift we derive the radio power 
at 1.4 GHz, Pi. 4, which we plot as a function of galaxy 
absolute magnitude in the ^'^r-band in the left-hand side 
panel of Fig. 2. Here, the dashed lines mark the area typ- 
ically populated by FRII-type objects according to Owen 
& White (1991). The FRI/FRII classification was first in- 
troduced by Fanaroff & Riley (1974): in FRI sources the 
radio emission peaks near the center of the galaxy and the 
emission from the jets fades with galacto-centric distance. 
FRII sources, on the other hand, have edge-brightened ra- 
dio lobes. As is evident from Fig. [2] the population of ra- 
dio galaxies in the local Universe is clearly dominated by 
FRI-type objects, which are typically hosted by elliptical 
galaxies with very little on-going star formation and weak 
emission lines (Ledlow & Owen 1995; Govoni et al. 2000; 
Best et al. 2005b). In the right-hand panel of Fig. 2 we 
plot Pi. 4 of our sample galaxies as a function of their stel- 
lar mass A/,; the solid line indicates the locus where the 
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Table 2. Optical Activity of Radio Galaxies. 





SF 


CO 


AGN 


NOA 


LRA centrals 


685 


606 


653 


651 


LRA satellites 


149 


91 


79 


71 


HRA centrals 


322 


426 


538 


1308 


HRA satellites 


66 


61 


57 


182 



'^/ 



Notes: The numbers of LRA and HRA centrals and satellites that 
are also identified as star-forming galaxies (SF), composite galax- 
ies (CO), optical AGN, or that are not optically active (NOA). 



specific radio power, Pn/M, = 10^^^^WHz"^Mq\ which 
coincides with the peak of the corresponding distribution 
(cf. Fig. 15). We split our radio galaxies in sources with 
low radio activity (LRA) and those with high radio activity 
(HRA), depending on whether P1.4/M* is smaller or larger 
than 10^^/i^WHz"^Mq\ respectively. The numbers of LRA 
and HRA galaxies identified in our sample, and split accord- 
ing to centrals and satellites, are indicated in Table 1. 

The optical and radio activity defined here are not mu- 
tually exclusive. In fact, it is well known that radio galaxies 
can also display optical activity. This overlap is summarized 
in Table 2, where we list the number of LRA and HRA cen- 
trals and satellites with optical activity (SF, CO or AGN) 
and with no optical activity (NOA). Note that the majority 
(75 percent) of the LRA galaxies also reveal optical activity, 
roughly equally split among SF, CO and AGN. However, in 
the case of the HRA galaxies 50 percent does not reveal any 
optical activity. Forty percent of HRA galaxies with optical 
activity are identified as AGN, while only 26 percent show 
evidence for ongoing star formation. 



3 THE ECOLOGY OF GALAXY ACTIVITY 

Using the sample defined above, we now investigate how the 
optical and radio activity of galaxies correlate with their 
stellar mass and with the mass of the dark matter halo in 
which they reside. In addition, we also examine whether the 
group hierarchy (central vs. satellite) has any impact on the 
presence of activity in the optical and/or the radio. 

We define f{A\P) as the fraction of galaxies with prop- 
erties P (halo mass, stellar mass, group hierarchy) that are 
in activity class A (SF, CO, AGN, LRA, HRA). These frac- 
tions are computed using 



Np 



f[A\p)= Y.'"^/Y.' 



(1) 



Here Wt = 1/T4nax,i is the weight of galaxy i, defined as the 
reciprocal of the comoving volume of the Universe out to 
a comoving distance at which the galaxy would still have 
made the selection criteria of our sample, Na\p is the num- 
ber of galaxies with properties P that are in activity class A, 
and Np is the total number of galaxies with properties P. 
Errors are determined using the jackknife technique. We di- 
vide the group catalogue into iV = 20 subsamples of roughly 
equal size, and recalculate f{A\P) 20 times, each time leav- 
ing out one of the 20 subsamples. The jackknife estimate of 
the standard deviation then follows from 



\ 



JV-1 

N 



E(/-/T 



(2) 



with fi the fraction obtained from jackknife sample i, and 
/ the average. 



3.1 Halo Mass Dependence 

The upper and lower panels of Fig. [3] show /(^|Mh,C) and 
/(j4|Mh,S), respectively. Here 'C and 'S' refer to centrals 
and satellites, respectively, A is the activity class (SF, CO, 
AGN, LRA, HRA), and Mh is the halo mass. Thus, these 
plots show the halo mass distributions of the various activ- 
ity classes, split according to centrals and satellites. Note 
how the peak of the halo mass distribution of centrals shifts 
towards a higher halo mass as their activity changes from 
star formation to optical AGN activity to radio emission. 
Star-forming centrals prefer the least massive environments 
with a halo mass of ~ 10^^'^h~^ Mq (this is the minimum 
halo mass probed by the galaxy group catalogue of Y07); 
their fraction drops rapidly with increasing halo mass, and 
becomes negligible in haloes with M >, 10^*/i~^ Mq. Cen- 
tral galaxies classified as composites prefer somewhat more 
massive haloes; the halo mass distribution peaks around 
~ li)^^h~^ M0 and remains non-negligible up to the scale of 
massive clusters. The fraction of centrals harboring an opti- 
cally classified AGN rises sharply from Mh = 10"'®/i"^ M© 
to lO^^'^ft"^ M0, stays roughly constant at a level of about 
14% out to Mh ~ 10^^-*/i"^ M0 and decreases by nearly a 
factor of 2 at the massive end. Finally, the percentage of cen- 
tral radio galaxies increases with increasing halo mass, and 
reaches a maximum in the most massive halo mass bin. The 
results for HRA and LRA centrals are remarkably similar, 
with a weak trend for the HRA centrals to prefer more mas- 
sive haloes than LRA centrals. Thus, when moving to more 
massive haloes, the characteristic activity of central galax- 
ies changes from star formation to optical AGN activity to 
radio activity. Note, though, that the absolute fraction of 
'active' galaxies decreases quite drastically with increasing 
halo mass; while more than 60 percent of the centrals in low 
mass haloes (Mh ~ lQ^^'^h~^ M©) are active (mostly star 
forming), at the massive end only ~ 20 percent of the cen- 
trals are active, split roughly equally among optical AGN 
activity and radio activity (both HRA and LRA). 

With regard to satellite galaxies, the fraction of those 
with star-formation activity diminishes by a factor of 3 going 
from group masses of -- 10"''/i"^ M© to lO^'^/i"^ Mq. Note 
that in massive haloes, the fraction of satellite galaxies with 
star formation activity is an order of magnitude higher than 
that of central galaxies. The fractions of CO and AGN satel- 
lites increase with increasing halo mass. Especially for the 
composites, this behavior is quite different from that of the 
centrals. Finally, the fractions of HRA and LRA satellites is 
rather constant across the full range of halo masses probed. 
Note, though, that only a very small fraction of satellite 
galaxies ( ^ 0.2 percent) are radio active, which is roughly 
an order of magnitude lower than for central galaxies. 

As discussed in Section 12.21 optical and radio activity 
are not mutually exclusive (cf. Table 2). Panel (a) of Fig. [3] 
shows the distribution of Pi. 4 of all radio sources for each 
optical activity class. Here, the fractions on the y-axis are 
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Figure 3. The activity fractions /(j4|Mij,C) for centrals (upper panels) and /(j4|Afij,S) for satellite galaxies (lower panels), as functions 
of their halo mass, M^^. These fractions reflect the halo mass distributions of the various activity classes (SF, CO, AGN, LRA and HRA, 
as labeled). The errorbars have been computed using the jackknife method, as described in the text, and correspond to Icr. 



computed with respect to the total number of radio galax- 
ies in our sample. SF galaxies dominate at Pi. 4 ^ lO^^*^ 
WHz~^, while the AGN distribution extends to larger val- 
ues, in agreement with the findings of Best (2004). The 
NOA radio sources are, instead, preferentially associated 
with P1.4 ^ 10^^ WHz"\ Panels (b) and (c) show the frac- 
tions f(Ac,pt + ^radio|Afh, C) , where Aopt corresponds to an 
optical activity class (SF, CO, AGN or NOA) and Aradio to 
a radio activity class (HRA or LRA). Upon inspection, these 
fractions are qualitatively consistent with what one would 
expect if radio and optical activity are independent, i.e., 

/(Aopt + A,adio|Mh, C) = /(v4opt|Mh, C) X /(yl,adio|Mh, C). 

Thus, for example, central radio galaxies with star forma- 
tion are typically found in less massive haloes than cen- 
tral radio galaxies with optical AGN activity, irrespective of 
their specific radio power. In the case of satellite galaxies, 
the statistics are typically too poor to reach any meaningful 
conclusion. This finding that the optical and radio activi- 
ties of galaxies are independent is in good agreement with 
Best et al. (2005a), who found that the probability that a 
galaxy of given stellar mass classifies as HRA is independent 
of whether it it optically classified as an AGN. 



3.2 Stellar Mass Dependence 

We now turn our attention to the stellar mass dependence 
of galaxy activity. The upper and lower panels of Fig. [S] 
show /(A|Mt, C) and /(y4|M,, S), respectively, where M* is 
the stellar mass of the galaxy. Note how the typical stel- 



lar mass of active centrals increases going from star form- 
ing centrals to centrals with optical AGN activity to radio- 
active centrals. The stellar mass distribution of SF centrals 
peaks at ~ 10^'^h~^ M©, and drops rapidly with increasing 
M*, becoming negligible for M, > 10"/i"^ Mq. The stel- 
lar mass distributions of CO and AGN centrals peak near 
10^°'' /i"^M0, though the distribution of the AGN centrals is 
somewhat more skewed towards larger masses. Finally, cen- 
tral galaxies with radio activity are typically very massive, 
with stellar masses in excess of ~ W^^h~^MQ. Qualitatively, 
these trends are remarkably similar to those with halo mass 
shown in Fig. |31 This similarity owes to the fact that the 
stellar mass and halo mass of central galaxies are tightly re- 
lated (e.g., Yang, Mo & van den Bosch 2008), and indicates 
that it will be difficult to disentangle whether the activity 
of central galaxies is causally connected to halo mass or to 
stellar mass (see below). 

As for the satellites, the stellar mass distributions are 
remarkably similar to those of the centrals. This is very dif- 
ferent from the halo mass distributions (cf. Fig. |3]), which 
show clear differences for centrals and satellites. This is a 
first indication that the activity of a galaxy is more closely 
related to stellar mass than to halo mass. 

Finally, we note that the stellar mass distributions for 
various activity classes shown in Fig. [S] are in good agree- 
ment with previous studies (e.g., Brinchmann et al. 2004; 
Kauffmann et al. 2004; Best et al. 2005b; Kewley et al. 2006). 



8 A. Pasquali et al. 



12 


_i_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i_y 




_ a 


1 — 1 


_ 


10 


- 




ALL - 


6^ 


I 




I 


fi 


: 


^ 


n 


; 


•2 6 




-1_ 




— 


o 


- 


r^ 


— 1 


- 


£ 4 


— , 


— 


— 


— 1 


— 


&^ 


- 












- 


2 
n 


\ff 


4 


1 1 1 1 1 


— 


~Sz 


^1-1-1 1 T 



21 22 23 24 25 
Log.o(Pi.4 [W Hz-']) 



2 - 



1 - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1- 


- b 




Centrals: LRA 


- 


- 




- 


- 




- 




- 




- 


^li 




- 


-; Li-i-|-T^ 


1 1 1 1 ill 1 1 1 



1 - 



I I I I I I I I I I I I I I I I I I I I I 



HRA 




Q II m I I I I * I i-i-rH I I Q I I fcjj^ i I ~ I . I I . I 

11 12 13 14 15 11 12 13 14 15 

Log,„(M, [h-'Mj) Log,,(M, [h-'Mj) 



Figure 4. Panel a: the distribution in radio power of all sources, distinguished among star-forming galaxies (purple line), composite 
galaxies (blue line), AGN hosts (green line) and radio sources with no optical activity (black line). The fractions on the y-axis are 
computed with respect to the total number of radio sources in our sample. Panel h: the distribution of radio-emitting centrals with low 
radio activity (LRA) as function of halo mass, split among star-forming galaxies, composite galaxies, AGN hosts and sources with no 
optical activity (the colour coding is as in panel a) . The fractions on the y-axis are calculated with respect to the total number of central 
galaxies in our sample. Panel c: as in panel b, but for radio centrals with high radio activity (HRA). 
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Figure 5. The activity fractions /(A|A/»,C) for centrals (upper panels) and /(A|A/»,S) for satellite galaxies (lower panels), as functions 
of their stellar mass, M*. These fractions reflect the stellar mass distributions of the various activity classes (SF, CO, AGN, LRA and 
HRA, as labeled). 



3.3 Breaking the causality degeneracy 

In order to examine whether halo mass or stellar mass is 
more causally connected to galaxy activity, one needs to 
probe the halo mass dependence for fixed bins in stellar mass 
and vice versa. Fig. [6] plots the fractions /(A|Mh,M*,C) 
(left-hand panels) and f{A\Mh,Mt,S) (right-hand panels) 
as functions of halo mass and stellar mass. Only (Mh, M»)- 



bins with at least 50 galaxies are shown, with darker shades 
of grey indicating larger fractions. Note that this figure re- 
veals the same overall trends as in Figs. |3]and[5] which are 
projections of Fig.[6]along the M, and Mh axes, respectively. 

The relatively tight correlation between stellar mass 
and halo mass for centrals, which partly owes to the way 
that halo masses have been assigned to the groups (see Sec- 
tion [2] and Y07), limits the dynamic range over which we 
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Figure 6. The activity fractions /(yl|Mij, M, , C) (left-liand panels) and /(A|Mij, M», S) (right-hand panels) as functions of halo mass 
and stellar mass. Only bins with at least 50 galaxies are shown, with darker shades of grey indicating larger fractions. From top to 
bottom, results arc shown for the following activity classes: SF, CO, AGN, LRA and HRA, as indicated. For central galaxies, the stellar 
mass is tightly correlated with halo mass, preventing us from investigating which of these two parameters is more causally connected to 
the activity. For the satellite galaxies, however, the fractions are almost independent of halo mass at fixed stellar mass, indicating that 
stellar mass is more causally connected to galaxy activity than halo mass. 



can probe the halo mass dependence at fixed stellar mass, 
and vice versa. Hence, it is virtually impossible to tell which 
of these two conditionals is more causally connected to the 
activity. In the case of the satellites, however, it is clear that 
the fractions vary more with stellar mass at fixed halo mass 
than with halo mass at fixed stellar mass, irrespective of the 
activity class being considered. In fact, the fractions are al- 
most independent of halo mass at fixed stellar mass, which 



clearly indicates that stellar mass is a more defining prop- 
erty of the activity class to which the galaxy belongs than 
halo mass. 

To make this more quantitatively, we ignore the distinc- 
tion between centrals and satellites and define the reduced 
conditional fractions 



Qi 



f(A\M„M^) 
fiA\M,) 



(3) 
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Figure 7. The reduced conditional probability functions, Qi (left-hand panels) and Q2 (right-hand panels), as functions of both stellar 
mass and halo masses. Only bins with at least 50 galaxies are shown, with darker shades of grey indicating larger values of Q. No 
distinction is made between central and satellite galaxies. From top to bottom, results are shown for the following activity classes: SF, 
CO, AGN and RADIO (LRA and HRA combined), as indicated. By definition, any variation of Qi with M-^ at fixed M^ indicates a 
dependence of galaxy activity on halo mass. Conversely, a variation of Q2 with M« at fixed M^ indicates a dependence of galaxy activity 
on halo mass. It is clear from this plot that stellar mass is more tightly related to galaxy activity than halo mass. 



and 



Q2 



f{A\hh,M,) 
f{A\M^) 



(4) 



Note that any variation of Q\ as a function of Mh at fixed 
Mt traces the dependence of galaxy activity on halo mass; 
conversely, changes in Q2 as a function of M* at fixed Mh 
reflect that galaxy activity depends on stellar mass. The 
reduced fractions Qi and Q2 are plotted in the left-hand 
and right-hand panels of Fig. [T] respectively, as a function of 
halo mass and stellar mass for the different activity classes. 
Note that we have combined the HRA and LRA classed 
into a single class (labeled RADIO) in order to increase the 
statistical significance of the present analysis. As in Fig. |6l 
darker shades of grey correspond to larger values of Qi and 
Q2. 

While at a given stellar mass, Qi varies little with 
halo mass, Q2 changes dramatically with stellar mass at 
a given Mh. This confirms that galaxy activity is more 
causally connected to stellar mass than to halo mass. For 
the star forming galaxies and the optical AGN, we find 
that the Q2(log A/, | logMh) are well fitted by simple linear 
relations. The right-hand panel of Fig. |8] plots the corre- 



sponding slopes, dQ2/dlog A/*, as functions of A/h. Clearly, 
these differ significantly from zero, with little dependence 
on halo mass, indicating that the fraction of galaxies with 
SF and optical AGN activity depends strongly on stel- 
lar mass at fixed A/h- For comparison, the left-hand pan- 
els plots dQi/dlogA/h as functions of M^ for the same 
activity classes. Although significantly different from zero 
(dashed line), these slopes are much smaller. The average 
dQi/dlog Afh is -0.25 and -0.20 for the SF and AGN galax- 
ies, respectively, indicating that there is a weak trend for 
SF and AGN activity to decrease with increasing halo mass. 
On the other hand, the average values of dQ2/dlog Af« are 
—0.80 and 2.0 for SF and AGN galaxies, respectively, indi- 
cating that while SF activity decreases with increasing stel- 
lar mass, optical AGN activity is strongly boosted in more 
massive galaxies. Comparing these values, the stellar mass 
dependencies of SF and optical AGN activity are ~ 3 and 
^ 10 times more prominent than the corresponding halo 
mass dependencies. This is the quantitative confirmation 
that the occurrence of galaxy activity is more strongly re- 
lated to galaxy stellar mass than to halo mass. 
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Figure 9. The central fractions, fccn as functions of stellar mass 
for all gala.xies (black line), and for four activity classes, as in- 
dicated. Note that being active increases the probability for a 
galaxy to be a central, indicating that being located at the center 
of the halo boosts both the star formation and AGN activity. 



3.4 Group hierarchy 

Our analysis thus far has clearly shown that stellar mass 
is a more important parameter for determining the activity 
class of a certain galaxy than halo mass. We now investigate 
whether the group hierarchy plays a role; does the mere 
fact that a galaxy is a central galaxy, rather than a satellite 
galaxy, boost or suppress its probability to be active in the 
optical and/or radio? 

To address this question we have computed the central 
fractions 



iVcon + iVaat 



Un{A,M,) = J2^-/ Yl 



(5) 



galaxies of activity class A and stellar mass Af*. Once again, 
the errors on the central fractions are determined using the 
jackknife technique with 20 subsamples. The results for SF, 
CO, AGN and RADIO (HRA and LRA combined) are shown 
in Fig.[9]as indicated. For comparison, we also show the cen- 
tral fractions of all galaxies, independent of their activity 
(or absence thereof). Clearly, the central fractions of active 
galaxies are higher than those of all galaxies, at basically all 
stellar masses, and for all activity classes. This indicates that 
galaxies of a given stellar mass are more likely to be active 
when they are centrals than when they are satellites. This 
efltect is most pronounced for star formation, while optical 
AGN activity seems to be only boosted by a mild amount 
for central galaxies compared to satellites of the same stellar 
mass. 

This is in qualitative agreement with Best et al. (2007) , 
who found that brightest cluster galaxies (i.e. equivalent to 
our 'centrals') are more likely to host a radio-loud AGN than 
other galaxies of the same stellar mass, and with van den 
Bosch et al. (2008a), who have shown that central galaxies 
have bluer colors, on average, than satellite galaxies of the 
same stellar mass. Since satellite galaxies of a given stellar 
mass where central galaxies of the same stellar mass prior 
to being accreted by their host halo, this suggests that ac- 
tivity, whether in the optical or in the radio, is suppressed 
(or quenched) as soon as a galaxy becomes a satellite. 



3.5 Halo-centric Distance 

The galaxy group catalogue also allows us to check whether 
the position of a galaxy within its group may foster or ham- 
per a specific kind of activity. For this purpose, we compute 
the projected distance, i?proj, of each satellite galaxy from 
the luminosity-weighted center of its group normalized by 
the characteristic radius of the group, Riso (see Section [2]). 
We define the normalized fraction of satellites of a given 



activity class and per bin of Rn 



i?proj/-Ri80 as: 



with Nccn and A'sat the numbers of central and satellite 



/no 



./ (-^norm |^, o) 
,f (-Rnorm|S) 



(6) 
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Figure 10. The solid lines show the radial distributions of satellite galaxies of different activity classes, corrected for the radial distribution 
of the general population of satellites unconstrained in activity (black solid line). Results are shown for star forming satellites (SF), 
composite satellites (CO), satellites with optical AGN activity (AGN), and satellites with no optical activity (NOA). If active satellites 
follow the same radial distribution as their parent population, they would follow the red horizontal line at /norm = 1- Clearly, CO 
and AGN satellites are more centrally concentrated than their parent population, while SF satellites have a preference to populate the 
outskirts of their parent haloes The dashed lines indicate the radial distribution of satellite galaxies that have the same stellar masses 
as those in the respective activity classes. The fact that these are very similar to the solid lines indicates that activity is not causally 
connected to halo-centric distance; rather the trends shown are due to mass segregation. See text for a detailed discussion. 



where /(-Rnorm|A, S) is the weighted fraction of satellite 
galaxies of activity class A in a given bin of -Rnorm, and 
/(-Rnorm|5') Is the Weighted fraction of satellite galaxies (un- 
constrained by activity) in the same _Rnorm-bin. This dis- 
tribution is therefore normalized with respect to the radial 
distribution of the general population of satellite galaxies. 

The solid curves in Fig. [10] show /norm as function of 
^norm, thus obtained. Results are shown separately for star 
forming satellites (SF; upper left-hand panel), composite 
satellites (CO; upper right-hand panel), satellites with opti- 
cal AGN activity (AGN; lower left-hand panel), and satellite 
with no optical activity (NOA; lower right-hand panel). No 
results are shown for radio activity, simply because the num- 
ber of satellite galaxies with radio activity is too small to 
allow for a meaningful assessment of their radial distribu- 
tion. As before, the (la) errors have been computed using 
the jackknife method. 

The SF satellites have /norm < 1 (> 1) at small (large) 
iinorm- This suggests that star formation is suppressed at 
small halo-centric radii. CO and AGN activity, on the other 
hand, seems to be boosted at small halo-centric radii. Al- 
though these results seem to indicate that halo-centric ra- 
dius plays an important causal role for galaxy activity, it is 
important to realize that we have not used stellar mass as a 
control variable. As shown in van den Bosch et al. (2008b), 
the group catalogue of Y07 used here reveals clear evidence 
for mass segregation, in that more massive satellites pref- 



erentially reside at smaller halo-centric distances. Since the 
analysis in Section [3.31 indicates that there is a strong causal 
connection between stellar mass and galaxy activity, the 
trends with halo-centric radius shown in Fig. llUl do not nec- 
essarily imply that activity is causally connected to -Rnorm- 
In order to test this we proceed as follows. For each satel- 
lite galaxy of activity class A in a given bin of i?norm we 
select a random satellite, unconstrained in activity, in the 
same radial bin (but not necessarily in the same group) that 
is matched in stellar mass (AlogM» $; 0.05) and redshift 
{Az ^ 0.01). We use this matched sample to compute 



f' 

J norm 



/a-/ (-Rnorm I S) 
/ (-f^norm |^^ j 



(7) 



where /M(finorm|S) is the fraction of satellite galaxies 
(unconstrained in activity) with halo-centric radius -Rnorm 
that are matched in stellar mass and redshift to those 
in /(-Rnorml^, S). The results thus obtained are shown in 
Fig.[lQ]as dashed lines. Overall, it is clear that /norm(^norm) 
is very similar to /norm(-Rnorm), for all activity classes shown. 
This indicates that the trends with halo-centric radius of the 
latter are simply a reflection of the mass segregation com- 
bined with a causal relation between activity and stellar 
mass. In other words, there is no indication that the oc- 
currence of galaxy activity in satellite galaxies has a direct 
dependence on halo-centric radius. 
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Figure 11. Left-hand panels: the distributions of log [L (Ha) /A/*] for star-forming centrals and satellites, normalized by their total 
number. Middle panels: scatter plots of the specific Ha luminosity as a function of stellar mass for centrals (upper panel) and satellites 
(lower panel). The solid red circles indicate the mean log[L(Ha)/Af*] per bin of stellar mass, while the errorbar indicates the error on 
the mean. Right-hand panels: same as in the middle panels, but as function of halo mass, Mjj. 



4 THE ECOLOGY OF GALAXY ACTIVITY 
STRENGTH 

In the previous section we have shown that the occurrence 
of any kind of galaxy activity is mainly governed by stellar 
mass, rather than halo mass. In addition, there is an indi- 
cation that being a central galaxy, rather than a satellite 
of the same stellar mass, causes a significant boost in the 
occurrence of activity. 

We now investigate the dependence of the activity 
strength on stellar mass, halo mass and group hierarchy. For 
star formation and optical AGN activity, we define the activ- 
ity strength in terms of the specific luminosities in the Hq 
and [QUI] emission lines, L(Hq:)/M* and L[OIII]/M*, re- 
spectively, both corrected for dust obscuration as described 
in Section 2. Since the spectra are obtained with fibers with 
a diameter of 3", we need to correct L(IIq) for the missed 
fiux. To that extent we scale //(Ha) by the ratio measured 
in the r band between the galaxy Petrosian flux and the flux 
observed through the spectroscopic flber. Note that no such 
correction is needed for I/[OIII], since the [QUI] emission is 
believed to be dominated by emission from the AGN. 

The results for the SF strength I/(HQ:)/Mt are shown 
in Fig. 111! The histograms in the left-hand panels show the 
distributions of log[Z/(Ha)/M«] for the SF centrals (upper 
left-hand panel) and SF satellites (lower left-hand panel), 
normalized by the total numbers of centrals and satellites 
in our group sample, respectively. Both distributions peak 
around L(Ra)/M, ~ 10~'^ Lq /{h~'^ M©) and are weakly 
skewed towards smaller specific luminosities. No clear dif- 
ferences are evident between the distributions for centrals 



and satellites. The middle and right-hand panels of Fig. [Tl] 
show scatter plots of the specific Ha luminosities as func- 
tions of stellar mass and halo mass, respectively. The solid 
red circles indicate the averages of log[I/(IIa)/M,] per bin 
of stellar mass or halo mass, while the errorbars (typically 
smaller than the solid circles) indicate the error on the mean. 
For central galaxies, the mean of log[L(Ha)/M,] decreases 
by a factor of four with increasing stellar mass and increas- 
ing halo mass. These trends are statistically significant at ~ 
4(T. Thus, in addition to suppressing the occurrence of SF 
activity, a larger stellar mass or a larger host halo mass also 
suppresses the strength of the SF activity. Unfortunately, the 
sample is not large enough to examine whether halo mass 
or stellar mass are more causally linked to this trend. For 
satellite galaxies, we find no obvious trend of the specific Ha 
luminosity with halo mass or stellar mass. This difference 
between centrals and satellites suggests that group hierar- 
chy has a direct impact on the strength of star formation 
activity in galaxies. 

Fig. [12] shows the same as Fig. [Tl] but for the specific 
[OIII] luminosity (related to optical AGN activity) rather 
than the specific Ha luminosity. Note that the distributions 
of log (1/ [OIII] /M,) are strongly skewed towards higher val- 
ues, for both centrals and satellites. In fact, the distributions 
are well fit by a sum of two Gaussians, with their respective 
peaks at log (L [OIII] /Af,) ~ —5 and —3. We identify these 
two 'populations' with LINERS and type II Seyfert galax- 
ies, respectively (cf. Kewley et al. 2006). In order to probe 
these two populations separately, we (somewhat arbitrar- 
ily) split our optical AGN in LINERS and type II Seyfert 
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Figure 12. Same as Fig. 1111 but for tlie specific [Ollf] luminosity of optical AGN. The vertical (horizontal) red line in the left-hand 
(middle and right-hand) panels corresponds to log(L[OIII]/M*) = —3.8, and indicates our split in LINERS (log(L[OIII]/-A/*) ^ —3.8) 
and type II Seyferts (log(L[OIII]/Af,) > -3.8). 



galaxies depending on whether log(I/[OIII]/A/*) is smaller 
or larger than —3.8, respectively. As is evident from the up- 
per panels in the middle and right-hand columns of Fig. 1121 
the central galaxies with Seyfert II activity reveal a decrease 
of their average specific [OIII] luminosity by a factor three 
going from low mass (either stellar or halo) to high mass 
(at the 4(T significance level). For the central LINERS, the 
(average) specific [OIII] luminosity diminishes by a factor of 
about 2 for increasing stellar and halo mass (at the 6(7 sig- 
nificance level). On the contrary, satellite galaxies show no 
dependence of optical AGN strength on either stellar mass 
or halo mass; neither for type II Seyferts nor for LINERS. 
Finally, Fig. [T3] shows the same as Figs. [TT] and [TJ] but for 
the specific radio power P1.4/M, [in W Hz'V^i'^M©)]. As 
evident from the left-hand panels, the distributions of spe- 
cific radio power for both centrals and satellites peak at 
log(Pi.4/Mt) = 12 (indicated by a red, vertical line), and we 
have used this value to split the population of radio galaxies 
in high radio activity (HRA; log[Pi.4/M*] > 12) and low 
radio activity (LRA: log[Pi.4/Af,] 5C 12) subsamples. The 
mean specific radio power of HRA centrals increases mildly 
with stellar mass by a factor of 1.5 (at the 3(t significance 
level), while the mean specific radio power of LRA centrals is 
constant with M* . In the case of radio satellites (both HRA 
and LRA), the mean log(Pi.4/M,), although noisy, shows 
no clear dependence on stellar mass. Similar trends are seen 
as function of halo mass (right-hand panels of Fig. llSfl . 



5 CONCLUSIONS 

An important aspect of galaxy formation is understanding 
the processes that trigger, quench and control star formation 
and AGN activity. Since both require the presence of gas, 
this is similar to understanding how galaxies accrete their 
gas, and how gas can be expelled, stripped or exhausted. A 
large number of processes are thought to play a role. At the 
low mass end, it is believed that supernova feedback plays a 
crucial role in expelling and/or heating gas so as to keep the 
overall star formation efficiency low. At the massive end, it is 
generally believed that AGN activity somehow expels or re- 
heats cold gas causing a quenching of star formation. Ram 
pressure and tidal forces are believed to remove gas from 
satellite galaxies, while major mergers are efficient in fun- 
neling gas to the central regions, igniting starbursts and/or 
AGN activity, which in turn may exhaust and/or expel the 
gas resulting in a subsequent shutdown of galaxy activity. 
In order to probe the environments in which these various 
processes operate, we have used the SDSS galaxy group cat- 
alogue of Yang et al. (2007) and examined the dependence 
of galaxy activity on stellar mass, halo mass, and group hi- 
erarchy (centrals vs. satellites). In particular, we have split 
the galaxy population in star formers, optical AGN, compos- 
ite galaxies (both star formation and optical AGN activity) 
and radio sources (split in high and low radio activity) us- 
ing emission line fiuxes extracted from the SDSS spectra by 
Kauffmann et al. (2003b) and by cross correlating the SDSS 
with the FIRST and NVSS catalogues. 

For the subsample of central galaxies, we find a smooth 
transition in halo mass as the activity of central galax- 
ies changes from star formation to optical AGN activity 
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Figure 13. Same as Fig. Illl but for the specific radio power of radio galaxies. The vertical (horizontal) red line in the left-hand (middle 
and right-hand) panels corresponds to log(Pi.4/M,) = 12, and indicates our split in high radio activity (HRA; log(L[OIII]/Af») > 12) 
and low radio activity (LRA; log(L[OIII]/A/») ^ 12) galaxies. 



to radio emission (see Fig. [Sj. Star-forming centrals pref- 
erentially reside in haloes less massive than W^^h~^ M©, 
central galaxies with optical AGN activity inhabit haloes 
with a characteristic mass of ~ 10^^ h~^ Mq and centrals 
emitting in the radio preferentially reside in haloes with 
Mh >, 10^*/i~^ M0. A similar trend is present as function 
of stellar mass. Star formation is mainly occurring in cen- 
tral galaxies with M, <, 10^^ h~^ M©, while the hosts of 
optical-AGN and radio sources have typical stellar masses 
of lO^°-*/i-^M0 and IO^-^/i-^Mq, respectively (see Fig.fS]). 
Unfortunately, since the relation between stellar mass and 
halo mass for central galaxies is relatively narrow, we cannot 
discriminate whether this trend is causally linked to stellar 
mass or to halo mass. 

In addition to the occurrence of activity, we have also 
investigated the halo and stellar mass dependencies of the 
activity strengths. The activity strength of star forming and 
optical AGN centrals is seen to systematically decrease with 
both stellar and halo mass by a factor of 4 - 2 across the 
range 10^-* < M. < lO^-^h"^ Mq and 10^^ < Mh < 
10^* h~^ Mq, while the specific radio power of HRA cen- 
trals increases by a factor of 1.5. Again, the tight relation 
between A/* and A/h prevents us from investigating which 
is more causally connected to the activity strength. 

In hierarchical models of structure formation, dark mat- 
ter haloes and galaxies continuously grow in mass. We may 
therefore interpret the mass dependence of central galax- 
ies as a trend with time: central galaxies start out form- 
ing stars, but once their stellar mass exceeds a few times 
10^''^"^ Mq (or their halo mass exceeds ~ lO^^/i"^ M©), 
star formation is quenched, and optical AGN activity be- 



comes important. By the time the halo mass starts to ex- 
ceed a few time 10^^h~^ M© (or the stellar mass exceeds 
~ 10^^ h~^ Mq), the radio- mode of AGN activity starts to 
become important, while optical AGN activity is somewhat 
suppressed. Interestingly, the mass scale of this first transi- 
tion coincides with the mass scale where (i) the mass-to-light 
ratios of dark matter haloes are minimal (Yang, Mo & van 
den Bosch 2003; van den Bosch et al. 2007), (ii) dark mat- 
ter haloes start to develop an accretion shock at their virial 
radius, giving rise to a hot corona (Birnboim & Dekel 2003; 
Keres et al 2005; Birnboim, Dekel & Neistein 2007), and (iii) 
galaxies undergo a transition from being disk-dominated to 
being bulge-dominated (Kauffmann et al. 2003a). 

The following picture emerges: central galaxies in haloes 
with A/h ^ W^^h~^ M© accrete cold gas, which they ac- 
tively turn into stars. During this period supernova feed- 
back is believed to play an important role (Dekel & Silk 
1986) and to be responsible for an increase of the baryonic 
mass fraction that is converted into stars with increasing 
halo mass. If a galaxy undergoes a major merger during 
this epoch, it regrows a disk around its newly formed bulge 
(e.g., Kauffmann, White & Guiderdoni 1993). The tight re- 
lation between bulge mass and black hole mass suggests that 
this bulge formation is somehow coupled to the formation 
of a supermassive black hole at its center (Gebhardt et al. 
2000; Ferrarese & Merritt 2000). Once the halo mass exceeds 
~ 10^^ h~^ Mq, it starts to develop a rapidly expanding ac- 
cretion shock, which heats up any newly accreted gas to the 
halo's virial temperature. Somehow most of this gas is pre- 
vented from cooling, resulting in a quenching of the star for- 
mation of the central galaxy (e.g., Dekel & Birnboim 2006; 
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Cattaneo et al. 2006; Birnboim et al. 2007). From this point 
onwards, the galaxy cannot regrow a new disk after it under- 
goes a major merger, which explains the transition towards 
more bulge dominated systems. These mergers are also likely 
to drive any remaining gas to the center of the merger rem- 
nant, where it gives rise to a starburst (exhaust the cold gas 
supply) and to AGN activity (e.g. Mihos & Hernquist 1996; 
Sanders & Mirabel 1996). The latter may play an important 
role in reheating or expelling gas, and thus in quenching the 
star formation (e.g., Menci et al. 2005; Springel, Di Matteo 
& Hernquist 2005; Hopkins et al. 2006). Once the halo has 
grown a hot corona that is massive and dense enough (i.e., 
when Mh starts to exceed a few times lO^^/i~^M0), it starts 
to dissipate the energy output from the AGN, giving rise to 
radio lobes, i.e. the hot gas serves as a 'working surface' for 
the jet (e.g. Kauffmann, Heckman & Best 2008). This in 
turn gives rise to radio-mode AGN feedback, which is be- 
lieved to play an important role in preventing the hot gas 
from cooling and forming stars (e.g.. Bower et al. 2006; Cro- 
ton et al. 2006; Kang, Jing & Silk 2006; Monaco, Fontanot 
& Taffoni 2007). 

As to satellite galaxies, our results support a picture 
in which satellite galaxies, in general, have their activity 
quenched. This is supported by the fact that the fraction 
of active satellites is always smaller than the fraction of ac- 
tive centrals of the same stellar mass (see Figs. (5)1, and by 
the fact that the probability for an active galaxy to be a 
central galaxy is higher than that for random galaxy of the 
same stellar mass (see Fig.[9l. The fact that satellite galaxies 
only reveal a weak dependence of activity on halo mass (with 
more massive haloes showing somewhat smaller fractions of 
active satellites), indicates that the efficiency of this activity- 
shutdown mechanism must operate with roughly equal effi- 
ciency in haloes of all masses. This is in good agreement with 
van den Bosch et al. (2008a,b) who, using the same galaxy 
group catalogue as that used here, found that satellite galax- 
ies are redder and more concentrated than central galaxies 
of the same stellar mass, but that the magnitude of this ef- 
fect is virtually independent of halo mass. As discussed in 
van den Bosch et al. (2008a), this seems to favor strangula- 
tion (i.e., the removal of the satellite's hot gas reservoir) as 
the main quenching mechanism, rather than processes such 
as ram-pressure stripping which mainly operate in the dense 
ICM of massive clusters. 

The activity of satellite galaxies shows a strong depen- 
dence on stellar mass, very similar to that of central galaxies. 
This is consistent with a picture in which satellite galaxies 
were central galaxies prior to being accreted, thus 'inher- 
iting' their stellar mass dependence. We also find that the 
activity of satellite galaxies is correlated with their projected 
group-centric distance: while star forming satellites prefer- 
entially reside in the outskirts, satellite galaxies with optical 
AGN activity have smaller projected distances to their group 
centers than the average satellite. However, these trends do 
not reflect a causal connection between activity and halo- 
centric distance. Rather, they are simply a reflection of 
mass segregation (more massive satellites have smaller halo- 
centric distances) combined with the relation between activ- 
ity and stellar mass. The absence of a direct link between 
activity and halo-centric distance indicates that the activity- 
quenching mechanism(s) must operate on time scales longer 
than, or comparable to the dynamical time in the haloes. 



This is again consistent with the idea that the main quench- 
ing mechanism is strangulation, which seems to operate on 
a timescale of several Gyrs (e.g. McCarthy et al. 2008; Kang 
& van den Bosch 2008), and is also supported by the lack 
of correlation between the activity strength of satellites and 
either stellar or halo mass as shown in Figs. 11, 12 and 13. 
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